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Abstract 



We introduce a mathematical setup for charge transport in quantum pumps 
connected to a number of external leads. It is proved that under rather general 
assumption on the Hamiltonian describing the system, in the adiabatic limit, the 
current through the pump is given by a formula of Biittiker, Pretre, and Thomas, 
relating it to the frozen 5-matrix and its time derivative. 

1 Introduction 

Transport in quantum pumps has been investigated in relation to various properties and 
from many perspectives |8| iH [H 123 113 CH] • The goal of this article is to provide a rigorous 
setting for a single but important aspect of these devices, namely the charge transport 
or, more precisely, its expectation value. The idealized setting is as follows: a pump, 
whose internal configuration varies slowly in time in a prescribed manner, is connected 
to n leads, or channels, along each of which independent electrons can enter or leave 
the pump. We assume that the electron in the lead has no transverse or spin degrees of 
freedom and may be thought of as a (non-relativistic) particle moving on a half line^. The 
incoming electron distribution, at zero temperature, is a Fermi sea with Fermi energy /i 
common to all leads. As a rule, this does not apply to the distribution of the outgoing 
electrons, as their energies may have been shifted while scattering at the pump. Because 
of this imbalance a net current is flowing in the leads. The expected charge transport is 
expressed by the formula [3 EI 



^Such extra degrees of freedom can be represented by adding channels. In general, the different 
channels may then have different propagation speeds. 





Here S = (Sij) is the n x n scattering matrix at energy /i computed as if the pump were 
frozen into its instantaneous configuration. A change of the configuration is accompanied 
by a change S ^ S + dS of the scattering matrix and by a net charge dQj leaving the 
pump through lead j. Finally e is the electron charge, which is henceforth set equal to 1. 

We shall next present a mathematical framework in which (jl.lj) can be phrased as a 
theorem. 

• The single-particle Hilbert space is given as 

n = no®L\R+X'"), (1-2) 

where states in L^(]R+,C") = ©"^]^L^(R_|_), resp. in Ho, describe an electron in one 
of the leads j = 1, . . . .n, resp. in the pump proper. The latter Hilbert space is not 
further specified, but the hypothesis A2 below confers on the pump the role of an 
abstract finite box [23]. Let Uj : H ^ H denote the projection onto Hq for j = and 
on j-th copy of L^(R+) for j = 1. . . . n. 

• Since the pump configuration is supposed to change slowly in time t, we will eventually 
consider the evolution of the electrons in an adiabatic limit, where s = et is kept fixed 
as e > tends to 0. In terms of the rescaled time coordinate s, called epoch, the 
propagator Us{s, s') on H satisfies the non-autonomous Schrodinger equation 

idsU,{s,s')=e-'H{s)U,{s,s') , (1.3) 

where H{s) is a family of self-adjoint Hamiltonians on Ti enjoying the following prop- 
erties: 

(Al) H{s) — H{s') is bounded and smooth in s , 

(A2) \\{H{s) + i)"™Ho||i < C for all s and some m e N 

(A3) H{s)iP = -d^tp/dx^ for ^ e C;^{R+, C") , 

(A4) app(i/(s))n(O,oo) = 0, 

(A5) H{s) = iJ_ for s < . 

Here || ■ ||i denotes the trace class norm over Ti, while the operator norm will be written 
as II ■ ||. Assumption A3 states that a particle in the leads is free; in particular, together 
with Al, it implies that changes in the Hamiltonian are confined to the pump proper: 

H{s)~H{s') = {H{s)-H{s'))Uo. (1.4) 

Assumption A4 requires that there are no positive embedded eigenvalues; A5 states 
that the pump is at rest for s < 0. 

• The initial state of the electrons should be an equilibrium state. This is achieved 
thanks to Assumption A5 by positing that the 1-particle density matrix at some (and 
hence any) epoch s_ < is of the form p{H_), where p(A) is a function of bounded 
variation with suppdp C (0, oo). A good example is the Fermi sea, where p(A) = 
6{E — A). The time evolution then acts as 

p{H.) ^ Ue{s, s.)p{H.)Ue{s., s) . (1.5) 
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• We define a generator of exterior scaling w.r.t. (jl.2j) by 

wliere v{x) : [0, oo) M is smootli witli v{x) = 0, resp. = x for small, resp. large 
X, and v'{x) > everywhere. We note that A = A* commutes with Uj, and set 
Aj = AUj. The operator A distinguishes between incoming and outgoing states, 
respectively associated with spectral subspaces A < —a and A > a with some large 
a > 0. Detection of a particle, and hence of its charge, deep inside lead j may be 
realized as the operator Qj{a) = f{Aj — a) + f{—Aj — a), where / G C°°(]R) is a 
switch function: /(a) = for a < —1, resp. = 1 for a > 1. The current operator then 
consistently is 

J,(a) = i[H{s), /(A, - a) + /(-A, - a)] =: /,+(a) + /,_(a) . (1.7) 

(In what follows we shall sometimes suppress the index j for the sake of notational 
simplicity). One feature of this choice of current operator is that the "ammeter" 
is located not at a fixed distance from the pump, but rather at a fixed number of 
wavelengths, a, from it: The longer the wavelength the more distant the "ammeter". 
In the case that one focuses on a narrow energy interval, say near the Fermi energy, 
the "ammeter" is also at essentially fixed distance from the pump. 

We remark that by the support property of v and by A3, the above commutator does 
not depend on s. The expectation value of the current at epoch s, i.e., in the state 
(jl.5|) . is then given as 

{I},{s,a,e) = tT{U,{s,s_)p{H_)U,{s_,s)I,{a)) . (1.8) 

In contrast to Qj{a), which clearly has an infinite expectation value in that state, Ij{a) is 
inclined to have a finite one. Moreover, it should behave as e if, in accordance with (ll.lj) . 
the charge dQj transferred during ds = e~^dt is to have a non trivial limit as e — > 0. 

Other realizations of the current operator are possible^, such as i[H{s), f{xj — a)], or 
the example based on the precession of a spin proposed in [TH] , and the result (jl.l|) should 
be independent of the choice. Our definition (jl.7|) has the property that Ij{a) naturally 
splits into two parts distinguished by their Heisenberg dynamics: The outgoing current 
Ij+{a) which is essentially free in the future and the incoming current Ij-{a) which is free 
in the past. As the initial condition is set in the past of the measurement, only Ij+{a) 
will be affected by scattering. 

• Finally we ought to state the reference dynamics Hq to which H{s) is compared in the 
definition of the (frozen) scattering operator [221 1211 

Sis) = SiHis), Ho) , (1.9) 

^The canonical choice of the current operators one normaUy finds in textbooks corresponds to f{x) = 

eix). 
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Figure 1: A scatterer with two channels (left) and its disconnected analog (right). Neu- 
mann boundary conditions are imposed at the edges of the channels. 

whose fibers S{s, E), called scattering matrices, appear in (jLip . While it turns out that 
the choice of Hq is largely irrelevant, for the sake of simplicity let Hq be the Laplacian, 
acting on L^(]R_,., C^), with a Neumann boundary condition at x = 0. S{s, k"^), {k > 0), 
then agrees with the familiar definition based on generalized eigenfunctions incident 
through lead i: {Sjie~^^^ + Sjie^^^)"'^^. In particular, this reproduces the standard form 



of the scattering matrix of two channels: 

S = 



r t' 
t r' 



with r, r' the right and left reflection amplitudes and t, t' the corresponding transmis- 
sion amplitudes. 

The fundamental equation may thus be given the following reformulation: 

Mm lime-'{I),is,a,6) = rdp{E)(^S*) . (1.10) 

In particular, for the Fermi sea p{HS) = 6{fi — HS) as initial state, we recover (ll.lj) from 
dp{E) = —6{E — p)dE. The limit a — >■ oo is taken so as to have the current measurement 
made well outside of the scattering region, but after the adiabatic limit e ^ 0. By doing 
so, a current is still measured within the same epoch as the scattering process which 
generated it, though at a different time. 

The result can be proven essentially in this form, though some problems that arise 
both in the infrared and the ultraviolet have to be dealt with. In fact, the adiabatic limit 
is realized jS] in the regime where the dimensionless quantity given by e times the dwelling 
time of an electron in the pump is small. Low energy particles have a large dwell time in 
the pumps and, in addition, may get trapped indefinitely as new bound states are born 
at the threshold E = 0. This means that no scattering description in terms of a single 
epoch is adequate at low energies. Similarly, at high energies resonances may become 
increasingly sharp with correspondingly long dwelling times ^. On the other hand, low, 
resp. high energy states do not contribute to the net current, since they are filled, resp. 



^If the pump is a chaotic bilHard, the classical dynamics will, in general, have arbitrarily long periodic 
orbits. It is natural to expect that such orbits give rise to resonances. 
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empty in both the incoming and the outgoing flow. We shall therefore concentrate on 
the contribution to the current coming from states in any intermediate energy range, as 
selected by the function x below. 

Theorem 1.1. Let x ^ C'o°(0, oo) with x = ^ on supp dp. Redefine the current operator 
(with a UV and IR cutoff) 

Ij,±{s,a) = x{H{s)y[H{s)J{±A,-a)]x{H{s)) , (a > 1) , (1.11) 

in 7| ). Then 

Inn lime-\l),{s,a,e) = rdp{E)(^S*) , (1.12) 

a-*oo e-*o Zn Jq \as / jj 

where S = S{s,E). The double limit is uniform in s ^ I, I being a compact interval, 
whence it carries over to the transfered charge /^^"^ dt'{I)j{et', a, e) = e^^ ds'{I)j{s', a, e). 

Remark. The auxiliary objects v and / affect the current operator. Nevertheless, they 
disappear from its expectation value in the adiabatic and the large a limits. This may be 
phrased as the statement that different ammeters measure the same current. 

We first give a heuristic derivation of Eq. ()1.12|) . which may however serve as a guide 
through the complete proof. (For more hints, see One could argue that for computing 
the expectation values of observables, which, like Jj-^_(s, a), pertain to the incoming part 
of phase space, scattering may be ignored, i.e., one may pretend that the state of the 
system at epoch s simply is p{H{s)). On the other hand, when discussing /j- _|_(s, a) one 
should do as if the state were 

piHis))+e[S^'\s),p{His))] + ... , (1.13) 

where 

/oo 
dtte'^^'^'H{s)e-'^^'^' . (1.14) 
-oo 

This comes from linearizing the Hamiltonian H{s + et) = H{s) + eH{s)t + . . . around 
epoch s, since U = 1 + eS^^\s) + ... is then the propagator from t = — oo to t = oo in the 
interaction picture based on H{s). Now (jl.l3|) is the state Up{H{s))U* after scattering. 
The contribution of the first term there to {I)j,+ cancels against {I)j,-, see Lemma 12.21 
The contribution of second term is responsible for the r.h.s. of Eq. ()1.12p . To foreshadow 
this outcome we formally rewrite it, see Lemma f5. 41 as 

[S«(s),p(i/(s))] = -id,,S{s',s),,=sP'{H{s)) , (1.15) 

with 

/oo 
c/te'^(^)*ij(s)e-'^(^)* . (1.16) 
-oo 

Here, 1 + (s' — s)dsrS{s' , s)s>=s + • • • is the Born approximation for the scattering ma- 
trix S{s',s) for the pair of autonomous Hamiltonians {H {s') , H [s)) . The importance of 
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Eq. (|1.15|) is twofold. First it reduces matters to "frozen" scattering data and their deriva- 
tives. Second, it makes it clear why only the variation dp{E) matters and, in particular, 
why for the Fermi sea p(A) = 9{fi — A) only states at the Fermi energy /i contribute to 
the current. 

Let us next discuss some of the conceptual issues involved in this result. 

1) Theorem 11.11 plays the role of the Kubo formula and like it gives a handle on 
transport, including dissipative transport, using Hamiltonian evolution. A well known 
weakness of linear response theory, stressed by van Kampen jljj , is that it takes the limit 
of linear response first and only then the thermodynamic limit. The correct order is, of 
course, the reverse. Theorem 11.11 is free of this criticism in that one starts with a system 
that has infinite extent, for which it is nevertheless possible to estimate the error made 
by linearizing the dynamics, see Lemma 14.11 The reason for that is the finite extent of 
the pump or, more precisely, the finite dwell times. 

2) Theorem 11.11 is an example of an adiabatic theorem for open, gapless systems. In 
contrast to other results of that kind, such as [HI El which establish that an embedded 
eigenstate evolves adiabatically, ours is about the evolution associated with an infinite 
dimensional spectral subspace, e.g. the Fermi sea. Moreover, the goal is not to establish 
that a spectral subspace is preserved by the dynamics; rather, it is to determine the 
amount by which it mingles with its complement, and the current is a measure thereof. 

3) On the l.h.s. of ()1.12|1 the scattering process and the measurement of the current 
are described as a single quantum history. This should be contrasted to usual textbook 
treatments of scattering rates, which have the following features: (i) The rates are com- 
puted classically as the product of a quantum mechanical scattering cross section and of 
an incident current; (ii) the use of the scattering cross section tacitly replaces the actual 
state by its free asymptote. For a static potential and a steady beam both steps have 
been justified |11^ in the limit where the measurement is deferred all the way to t —>■ oo. 
This is not quite so in the present fluctuating setting, where the ammeter is located at an 
intermediate scale, as specified by the order of the e, a limits: far away from the pump on 
the scale of the scatterer, but close to it on the scale of the distance traveled by a particle 
on the adiabatic time scale. 

4) The formula of Biittiker et al., Eq. ()1.1|1 . is pointwise in time. Suppose we represent 
the flow by means of coherent states of widths AE and At. In line with the remark just 
made, one might suspect that At needs to be picked small w.r.t. to the adiabatic time 
scale e~^, and hence AE ^ he hj the Heisenberg uncertainty relation. At temperature 

the 1-particle distribution of the incoming flow is (l + e^'-'^"^^)"^, which suggests that 
AE < f3~^. This seemingly presents an obstacle to applying scattering arguments when 
j3^^ < he. Our result, which applies to p(A) = 9{n — A), shows that this obstacle can be 
overcome, and that one can take the zero temperature limit (3 ^ oo before taking e — 0. 

Results related to Theorem 11.11 have been obtained at the level of both physical and 
mathematical rigor. Among the former we mention, besides of jHlEj, the formula of Lee, 
Lesovik and Levitov jl8|, which expresses the noise generated by quantum pumps (as 
well as other moments) in terms of the scattering data. However, the relations regarding 
noise are not local in time. Rather these are integral relations that hold for a cycle of 
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the pump. Physically, what sets the current apart from the noise is its linear dependence 
on p. Rigorous mathematical description of the adiabatic scattering of wave packets is 
discussed in [201 [1^]. Like our result, they rely on propagation estimates. 

The plan of the rest of the paper is as follows. In Sect.|21we shall verify that {I)j{s, a, e) 
and the Stieltjes integral in ()1.12|1 are well-defined. Further preliminaries, like propagation 
estimates, will be addressed in Sect. El In Sect. E]we shall compute the limit e — of the 
current in terms of data involving H{s) and H{s), but oblivious of the past {H{t)}t^s. The 
expression will further reduce to the "frozen" scattering data ()1.9|1 in Sect. El where the 
limit a ^ oo is taken. In the Appendix we establish some trace class estimates for opera- 
tors related to The main ideas are further discussed at the beginning of Sects. IHEl 

We conclude with a remark on notation. Multiple commutators are denoted by 
adJJ^^(5) = [aS^'^\B),A] with adf{B) = B. By F{A > a) we mean the spectral 
projection of A onto [a, oo). The trace class ideal is denoted as J'l. Generic constants are 
indicated by C. 

2 The current operator and the scattering matrix 

The state of independent quantum particles is described by a density matrix < p 
normalized so that tr p is the number of particles. In the case of Fermions, p < 1. 
Thermodynamic systems have tr p = oo and observables that are otherwise innocent, and 
bounded operators in particular, may fail to have finite expectation values. For example, 
the charge associated with the "box" f{A — a) in phase space is infinite. Nevertheless, the 
current flowing into the box should have a finite expectation. We begin by showing that 
the incoming and outgoing current operators of Eq. (|1.7p are trace class and consequently, 
the expectation value {I)j{s,a,e) in any Fermionic state is well defined. 

A classical interpretation of the trace class condition is for the observable to be asso- 
ciated with a localized (bounded) function in phase space. The heuristic reason why the 
current is trace class is then as follows: The commutator [H{s), f{Aj — a)] is localized near 
the boundary of the box: a curve (hyperbola) in phase space. This is where the ammeter 
is. The ultraviolet and infrared cutoff x{H{s)) then further delineate a compact region of 
phase space near the hyperbola. Our first preliminary result confirms this picture. 

Proposition 2.1. The operator of incoming and outgoing current in the j-th channel, 
Ijj-{s,a), is a trace class operator which is localized near ±a in the sense that 

\\F{\ATa\ > a)/,±(s,a)||i < C^(l + a)"^ (2.1) 

for all N en, a > 0. 

Proof. This proposition is a direct consequence of Lemma [3. 31 below. □ 

The pumping formula of Thm. 11.11 implies that no current flows if the pump is not 
operating, since dS/ds = if H{s) is independent of s. Indeed, in the state p{H{s)), like 
in a thermal state, different leads are at equilibrium with one another; moreover, in each 
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lead, right and left moving states yield compensating currents if equally occupied. When 
the pump is in operation, it will still prove useful to know that no persistent currents are 
flowing in the "adiabatic" state p{H{s)). 

Lemma 2.2. The currents in the state p{H{s)) vanish, namely 

tr{p{H{s)){Ij4s,a)+I,.{s,a)))=0. (2.2) 

Proof. We omit s from the notation. Since I±{a) G J7i it suffices to prove ()2.2j) for smooth 
p, since one may approximate the general case by a sequence s— lim^^oo Pn(-f^) = p{H). 

First we show that tr{p{H)I±{a)) is independent of a. Note that /±(ai) — /±(a2) = 
x{H)i[H,g{A,)]xiH) with g{Aj) = fi±A, - a,) - f{±A, - a^). Since x{H)g{Aj) = 
x{H)g{A)Iij e Ji, see fTOTj) . we have 

tr(p(i^)(/±(aO - /±(a2))) = Mpx{H)Hg{A,)xiH) - px{H)g{A,)Hx{H)) = (2.3) 

by cyclicity. It thus suffices to prove ()2.2|1 in the limit of a ^ oo. To this end, let Hq be the 
Neumann Hamiltonian on the leads introduced below Eq. ()1.9j] . and let J : L^(M_|_, C") —>■ 
Ti be the embedding given by p.2|) . We maintain that 

lim \\{p{H)I^{a) - Jp{Ho)liia)J*\\^ = , (2.4) 

a—i-oo 

where I±{a) is defined as in p. Ill) with H replaced by Hq. Indeed, by expanding that 
commutator we reduce matters to two estimates, both of the form 

lim \\xi{H)f{±A, - a)x2{H) - Jxi{H^)f{±A, - a)x2{H,)r\\^ = , 

a— »oo 

with Xi ^ C'^(]R), {i = 1,2). Then we write the difference as 

iXi{H) - Jxi{H,)r)f{±A, - a)x2{H) + Jxi{H,)r f{±A, - a){x2{H) - Jx2{H,)r) . 

Since Xi{H) - Jxi{Ho)J* G Ji by (EH) and s-lim„^oo /(±^j -a) = 0, Eq. (Q follows. 
At this point we only need 

tr(p(i/o)(4(a) + /°(a)) = 0, 

which holds by time reversal invariance: For KiIj{x) = tplx) we have KHqK = Hq, KAK 
= -A and hence KIl{a)K = -/° (a). □ 

Remark. One can establish the result without explicitly using time reversal and using 
instead the pull through formula Eq. (jA.lOj) . 

Mourre theory, see e.g. plays a double role in our analysis. First, it is at the heart 
of time dependent scattering theory and the propagation estimates that we shall discuss 
in the next section. At the same time, the theory also plays a role in time independent 
methods and we shall use it to establish the differentiability of the scattering matrix which 
appears on the right hand side of Eq. ()1.12|1 . 
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Proposition 2.3. Under the Assumptions A1-A5, the fibers of the frozen S matrix, 
S{s,E), are continuously differentiable in s at E > 0. In particular, the integral on 
the right hand side of Eq. il.l^) is well defined for p with bounded variation. 

Proof. The Hamiltonians H{s) are dilation analytic of type (A) with respect to the 
conjugate operator A. In particular, 

adf{H) is H{s)-hounded. (2.5) 

More importantly, for any energy E > the Mourre estimate holds 

Ea{H{s))i[H{s), A]Ea{H{s)) > eoEA{H{s)) (2.6) 

for some open interval A 3 E and > 0. Note that a compact term on the r.h.s 
can be dismissed by Assumption A4. The Mourre estimate ()2.6p is stable under small 
bounded perturbations of s, which is seen from ()A.5jl and from the fact that i[if(s),y4] 
is independent of s, just as the commutator in ()1.7j) . Therefore, ()2.(i|l holds uniformly in 
s G / and E E J , with J C (0, oo) compact, and so do the usual consequences of these 
assumptions. They comprise: 

i) Resolvent smoothness. Let {A) = (1 + A^)^ and r > 1/2. Then 

B{z,s) = {Ar'\H{s)-z)-\Ar 
has smooth boundary values aX z = E + iQ,{E E J)^ satisfying 

||a^5(E + iO,s)|| <Cfc (2.7) 

for A; = 0, 1, . . ., see [1^]. Since 

dsB{z, s) = -{A)-'iHis) - z)-'His)iHis) - ^)-^(A)-^ 

and H{s) = {A)^^H{s){A)~^ , the function B{E + lO, s) is jointly continuously dif- 
ferentiable in {E,s). 

ii) iY-smoothness. {A)^"^ is _f/'(s)-smooth as a consequence of ()2.7|) for k = 0. 

iii) Stationary representation of the scattering matrix. Let 

Ej{H{s))n ^ C^dE , ^ ^ {ij{E)} (2.8) 

be the spectral representation for H{s) on J, and set ro(-E)'?/' = ip{E). Then 
Vq{E){A)^''' : 7"^ — > C" is bounded by (ii). Let S{s',s) be the scattering opera- 
tor for the pair (H (s') , H (s)) . Its fibers S{E) : C" ^ C" with respect to ipiH) 
admit the representation 1^ 

SiE) = 1 - 27riro(A)(r - ViH{s') - (E + lO))" V)ro(A)* , (2.9) 

with V = H{s')—H{s). The r.h.s. is defined pointwise because of = {A)~'^V{A)~''' 
and of (i). From this, Eq. ()2.9|) . and the statements about B{E -\- lO, s) we see that 
S{s' , s)Ej{H{s)) is continuously differentiable in s' e /, and that the derivative 
can be computed fiber-wise. The same applies to the scattering operator S{s) = 
S{H{s),Ho), though ()2.9p appears slightly modified. 

□ 
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3 Propagation estimates 



Propagation estimates play a key and multiple role in our analysis. One role is that 
they guarantee that particles do not get stuck in the pump. Consequently, the scattered 
particle indeed sees a frozen scatterer to lowest order in the adiabatic limit, and a linearly 
changing scatterer to first order. At the same time propagation estimates also play a role 
in establishing that the current measured by the ammeter at epoch s is determined by 
the state of the pump at the same epoch, in the adiabatic limit. This may be interpreted 
as a statement that the particles neither linger nor disperse too badly in the channels 
Propagation estimates play yet another role in establishing the "rigidity" of the current: 
The expectation value for the current a-priori depends on a choice of a switch function 
/, the function v{x) in the generator of dilation and a choice of the initial configuration 
of the pump, H_. This dependence is suppressed from our notation because Thm. 11.11 
implies that the dependence disappears in the limit. Ultimately, this independence is a 
consequence of propagation estimates. 

Eq. ()2.(ij) implies a minimal escape velocity estimate jiT^ ITl] for the autonomous dy- 
namics generated by H{s). The constants involved are understood as being uniform in 
s G / and in the stated range for a. Within proofs we shall abbreviate H = H{s). 

Lemma 3.1. Let x in C^{0,oo) (in particular with support away from E = 0). Then, 
for some 6 > 0, for all a eR, b,t > and N e N: 

\\F{A <a-h + et)e-'^^'^'x{H{s))F{A > a)\\ < CN{e){h + Ot)-^ . (3.1) 

Similarly, if b,t < 0, then 

\\F{A >a-b + et)e-'^^'^'x{H{s))F{A < a)\\ < CNiO)\b + 9t\-^ . (3.2) 

Proof. The case of 6 = is covered by Thm. 1.1 in [14j, since its hypothesis (besides 
of (Eni)) that aS^\f{H)) is bounded for / G C^iM.) and A; > 1 holds true by (Q. To 
be precise, the result is formulated there for suppx C A, where A is as in (I2.(j|l and 
< 9 < 9q, but it extends to our case by a covering argument. Actually, the proof given 
there essentially covers the general case 5 7^ 0. More precisely, let 

Atr = T-\A - a + ^ - 9ot) 

and / G C°°(M) be a function with |/('=)(2;)| < C{x)-''J' < 0, and /(x) = for x > 0. 
Then 

\\f{A,^)e-^^\{H)F{A >a)\\< C^r'"^ 

uniformly in < t < r and a G M. For 6 = 0, this is equation (2.11) in ^3], whose proof 
applies to 6 > as well. Let r = h + t. Since 

h-9t 6/2 - 9ot , 

> - + 5 

h_+t - b + t 

"^Dispersion arises because we use a non-relativistic (quadratic) kinetic energy in the leads. In much of 
the hterature on pumps dispersion is circumvented by making the kinetic energy hnear in the momentum. 
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for some S > and all 6, t > we have 

for some / of the required type. □ 

The lemma after next is a refined version of Prop. 12.11 above. The discussion is sim- 
plified by the observation that the regularized currents Eq. (Il.ll|) may be written as 

I a) = xiHis))i[H,is), f{±A, - a)]xiHis)) , (3.3) 

where Hb{s) = H{s)b{H{s)) with b G C^(]R) and bx = X- The unregularized current, 
Eq. p.7p . is independent of s by A3 and the commutation [H{s),Ilj] = 0. Since the 
commutation fails when H{s) is replaced by -ffb(s) the regularized current is a-priori epoch 
dependent. However, the commutation is essentially recovered on the "box" f{±A — a) 
when a is large. More precisely: 

Lemma 3.2. For j = 0, . . .n and a > —1 we have 

\\[Hhis),Uj]f{±A-a){±A-a + i)''\\ <C (3.4) 

Proof. Since [A + i)^^/(±A — a){±A — a + i)^ is uniformly bounded in a > 1 we may 
prove instead that 

\\[Hbis),Uj]{A + if\\<C . (3.5) 
We begin by showing that IliHbIlj{A + i)'^ = IliHbIlj{Aj + i)'^ is bounded for i ^ j. Indeed, 

Hb{A, + if = {A, + ifHb + 2{A^ + i)[Hi„ A^] + [[i/,, A,], A,] , 

from which we infer that {Aj + i)~'^Hb{Aj + i)^ is bounded. Hence so is (use IliAj = 0) 



Il,HbIlj{A^ + lY = -U,{A, + i)-'Hb{A^ + lyu 



J ■ 



Now, setting H^- = 1 - H^- = XlLo.fcT^j ^^^^ 

[Hb, Uj]{A + 1)2 = {fl^HbU^ - U^HbU^YA + if , 
and the result follows. □ 

This implies that although the regularized current operator a-priori depends on the 
epoch and the pump, the dependence disappears in the limit of large a. 

Lemma 3.3. 

\\i[Hbis)Ji±A,-a)]xiHis))\U < C, (3.6) 

\\F{A<±a-a)i[Hb{s),f{±A,-a)]xiH{s))\\^ < 0^(1 + a)-'', (3.7) 

\\F{A>±a + a)i[Hf,is)Ji±A,-a)]xiH{s))\\i < C;v(l + a)"^ (3.8) 

for a, a > 1. The same bounds hold in operator norm if x{H{s)) is replaced with (±A — 
a + if. 
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The estimate (|3.8|) prevents the current operators from being located very far in the 
outgoing region of phase space. This will play a role in the next lemma and in Sect. 01 
Actually, instead of ()3.8|) _ we shall use there the weaker statement with characteristic 
function F{A > a + a). The pair and ()3.8|) _ keeps the current operators away from 

the pump, a property used in Sect. 

Proof. The estimates to be proven are of the form 

\\T[Ht„ f{±A, - aMH)\U < \\T[H,, f{±A, - a)]{±A -a + \\{±A -a + 

< C\\T[H,, f{±A^ - a)]{±A -a + if\\, (3.9) 

where we used fjA.16|) . 

We then have to establish the corresponding bounds for the remaining operator norm. 
Since g{Aj) = g (0) + {g (A) -g{0))U^ we have f{±Aj- a) = f{±A-a)Uj due to /(-a) = 
for a > 1. Writing / = //, where /(■) = /(■ + 2), the commutator in ()3.9|) is 

[H,, fUj] = [H,, f]U,f + fU, [H,, f] + f[H,, U,]f , 

with / = /(±A — a). In the contribution to (j3.9|l of the first two terms the projections 
Uj may be moved out to the right or to the left, using [T, 11^] = 0. As for the last term, 
we use fl3.4p . At this point, ()3.9p reduces to corresponding estimates for 

\\T[H,JiA)]iA + in + \\Tf{A)\\ , (3.10) 

where we replaced ±A — a by A, as the argument given below applies to the more general 
case. 



In the case of ()3.6p . where T = 1, the second term ()3.10p is clearly bounded. On the 
first term we use the commutator expansion, see eq. (B.16), or [10], based on the 
Helffer-Sjostrand representation of f{A). 

m—l ^ 

[H,J{A)] = ^-(-l)'=-iadJ)(/7,)/W(A) + i?^, (3.11) 

k=l 

Rm = j dxdyd-J{z){A-z)-^adf\H,){A-z)-^ , 

where = {d^ + idy)/2 and / is an almost analytic extension of /, which may be chosen 
in such a way that 

/m+2 
dxdy\d-J{z)\\y\-^-' < 
ifc=0 

for p = 0, . . .n, the norms being defined as ||/||fe = / dx{x)^\f{x)\. The choice of jT^ is 
such that if supp/' is compact, as it is in our case, then \y\ > Ci\x\ — C2, (Ci, C2 > 0), 
for z = X + ly & supp /, which thus implies 

\\{A-z)-\A + ,)\\<C{\y\-' + l) (3.12) 
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The expanded terms remain bounded if multiplied by (A + i)^ from the right. For the 
remainder we obtain 

m+2 

\\Rm{A + l?\\ < C\\e.df\H,)\\ Y,{\\f^'^\k-m-l+\\f''^\\k-m+l) , (3.13) 

k=0 

which is finite for m > 3. 

We shall next prove ()3.7|) with F{A < —a) replaced by F{A < —3a) for easing notion 
below. This may be further replaced by f{—a~^A — 2) because of F{A < —3a) = F{A < 
-3a)f{-a-^A - 2). Since /{-a'^A - 2)f{A) = 0, the claim Eq. (jTTnil reads 

\\[H,Ji-a-'A-2)]fiA){A + if\\<Cr,a-''. (3.14) 

We now apply the expansion ^ni} to -a-^A - 2 instead of A. Since f^''\-a-^A - 
2)f{A) = only the remainder contributes. To bound ||-Rm(^ + i)^ll II ("'^^^^ ~ 

2 — 2;)^^(A + i)|| < Ca{\y\^^ + 1) instead of ()3.12p and collect the powers of generated 
by each commutator. We so obtain the bound of ()3.13|) times q;~™+^. 

Finally, ()3.8j) can be brought into a form similar to ()3.7j) by replacing / with / — I. Both 
functions have the same commutator with Hb but essentially complementary supports. □ 

In Sect. E]we shall describe the dynamics ()1.3p in terms of the autonomous dynamics 
generated by the Hamiltonians H{s). Once this is achieved, the choice of an argument in 
initial 1-particle density matrix p, be it H_ or H{s), does not matter much for the current 
measurement, as defined by (jl.ll|) . This is the content of the following lemma. 

Lemma 3.4. Let p, x be as in Theorem \l.l\ Then, for a > 1, 

hm l|(p(ff(s)) -p(i7_))e-'^W*J±(s,a)l|i = . (3.15) 

t— >— oo 

We remark that no claim of uniformity w.r.t. s, a is made here. The statement for 
t +00 is also true, but not needed. 

Proof. We first consider the case where p is smooth, and since H is bounded below we 
may assume p G C^(]R). Then 

p{H) - p{H_) e Ji (3.16) 

by Eq. ()A.5|) . To estimate the trace norm in ()3.15|) we use ()3.3|) . insert 1 = F{A < 
a + a) + F{A > a + a) to the left of the commutator in ()1.11|) . and use ()3.6t 13. 8|) . We so 
obtain the bound 

II (piH) - p{H_))e-''''xiH)F{A <a + a)\\-C + C^v""^ . (3.17) 
The remaining operator norm is bounded for t < by 

\MH) - p{H^))F{A < a + a + 9t)\\ + C^r^ , (3.18) 
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where we used (13. 2|) with 6 = 0. We pick a = —6t/2, so that the remainder in (|3.17p 
tends to zero, and s— hm^^.oo . .) = in ()3.18|) . Since p{H) — p{H_) is compact, the 
norm vanishes in the hmit. 

In the general case, where p is of bounded variation, let p„ G C^(M), sup„ ||pn||oo < oo, 
be a sequence such that Pn(A) p{\) pointwise, whence s— lim„^oo Pn(-ff) = p{H), and 
the same for H_ instead of H. Also, s— limj^_oo(e~'^*Pac(-f^)~e~'^-*r2_) = by definition 
of the wave operator fl_ = Q_{H_, H) |22ll2ni- Since I±{s,a) is trace class by ()3.6|) and 
since 

s-limX„ = 0, YeJi^ lim ||X„r||i = , (3.19) 

n— >oo n— >oo 

we have 

hm \\{pn{H)-p{H))I±{s,a)\U = 0, 

n— >oo 

lim ||(p„(/7_)-p(/7_))^]_/±(s,a)||l = 0, 

n— >oo 

lim ||(e-'^* -e-'^-*fi_)/±(s,a)||i = 0, (3.20) 

t— »— oo 

where we dropped Pac(-f^) due to Psjc{H)x{H) = x{H). We then estimate {I± = I±{s, a)) 



'±||i 



||(p(i7)-p(iJ_))e-'^% 

-iHf.T / /■ TT\ ITT ,, 

'±l|l 



< \\{p{H) - p„(i7))e-^*/±||i + \\{p^[H) - p„(i/_))e-^*J, 



+ ||(p.(i7_) - p(i7_))(e-^^' - e-'^-'r]^)/±||i + ||(p„(ff_) - p(i^_))e-^^-*^]_J±l|l 
l|(p(/f) -p„(i/))/±||i + ||(p„,(/f) -p„(if_))e-^*/±||i 
+C||(e-^*-e-'^-*^]_)/±||l + ||(p„(iJ_)-p(iJ_))^]_/±||l 



Given e > we first pick n large enough such that the first and the last term together 
are less than e/2. For large negative t the same is true for the second and the third term 
by (|3.20|) and the first part of the proof. □ 



4 The adiabatic limit 

The current generated by adiabatic pumps can be understood within the general frame- 
work of the theory of linear response: The adiabatic limit corresponds to weak driving, 
and the formal derivation of Eq. (Il.lj) in [Hj is a perturbation expansion. Formally, the 
change in the state p of system, obtained by linearizing the Hamiltonian around epoch s 
is 

Uls,s.)p{E{s))Uls^,s)-p{E{s))^-^e T t e'^(^)*[if(s), p(i/(s))]e-^(^)*. (4.1) 

J — oo 

One immediate difficulty with this expression is that the integrand on the right hand side 
grows linearly with time. As an operator identity the above equation does not make sense, 
and the right hand side is not recognizably of order e. 

One of the reasons why a perturbation expansion can nevertheless be made is that the 
current of Eq. (jl.llj) has good localization in phase-space and so distinguishes a region 
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of phase space where p is to be evaluated. Therefore, only a restricted range of times 
contribute to the integral: The time associated with propagation from the pump to the 
ammeter. This makes the expectation value of the current well defined. Estimates of this 
kind are known as propagation estimates and are controlled by Mourre theory. 

In this section the limit e — in ()1.12j) is taken as a is kept fixed. The main result is 
the following. 

Proposition 4.1. For fixed a>l, we have 

hme-i tr([t/,(s, s.)p{H^)Ue{s^, s) - p{H{s))]I±{s, a)) = ii{P^l\s), p{H{s))]I±{s, a)) 

(4.2) 

uniformly in s & I , where 

Q^l\s) = -i dtte'^(^)*if(s)e-'^(^)* . (4.3) 

J —CO 

Moreover, 

lime-^ tr (f/e(s, s^)p{H^)Ue{s., s){I+{s, a) + /_(s, a))) = 

tr{[dl\s),p{Hismi+is,a) + I.is,a))) . (4.4) 

Remarks. 

• The integral ()4.3|) is a trace class norm convergent integral once multiplied from the 
left or from the right by I±{s, a), as in ()4.2j) . 



Eq. ()4.2|) separately describes, to leading order, the incoming and outgoing currents in 
an adiabatically evolved state as compared to the corresponding instantaneous state 
p[H{s)). Eq. ()4.4j) is then an immediate consequence of ()2.2j) . 

i^^'' (s) is formally the first order in e term in the expansion of the wave operator which 
for fixed s relates f/e(s', s) to the autonomous dynamics Us{s' — s) generated by the 
Hamiltonian H{s): 

u,{a) := e-'^^^)^^''^ , 

where e has been suppressed from the notation. In the next section we shall also meet 
the first order term in the expansion of the scattering operator p.l4j) relating these 
two dynamics. 



/oo 
citte'^W*if(s)e-'^(^)* 
-oo 



(4.5) 



The reason fi^''(s), rather than fi^''(s), appears in ()4.2|) is that the initial condition 
was set in the past of the current measurement. 

All estimates in this section hold true uniformly in s G /, s_ < 0. Constants may 
however depend on a. 
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The two dynamics, and Ug, are compared by means of the Duhamel formula 

Ue{s', S) = Us{s' -S)+ 16-^ j^^ dSiUeis', Si){H{Si) - H{s))us{si - s) . (4.6) 

Starting at epoch s, the Heisenberg dynamics of the currents carries them through the 
scatterer within a finite time, i.e., essentially still at same epoch s. When applying ()4.6|) 
to I±{s,a) it is thus appropriate to linearize H{si) — H{s) around s. The next lemma 
essentially says that calculating the current to first order in e the error one makes is 
second order in e: 

Lemma 4.2. We have 

e-'\\[U,{s.,s) - Ms. - s)+X,{s_,s))]I±{s,a)\\i < Ce , (4.7) 

where ^ 

Xe(s_, s) = ie~'^ / dsi{si - s)t4(s_, si)H{s)us{si - s) . 

J S- 

In turn this yields the following expression for the (rescaled) currents at epoch s: 

e-' tr(p(if_)[/,(s_, s)I^{s, a)U,{s, s.)) = 

e''^tr{us{s - s.)p{H^)us{s^ - s)I±{s,a)) + 
tr(X,(s_, s)*p(//„)m,(s_ - s)/±(s, a)) + 
e"^ tr{usis - s_)p(if_)Xe(s_, s)/±(s, a)) + 
e-^ tr(p(i7_)Xe(s_, s)/±(s, a)X,(s_, s)*) + 0e . (4.8) 

Proof. Using Duhamel's formula ()4.6|1 and 

H{si) - H{s) - (si - s)H{s) = {H{si) - H{s) - {s, - s)H{s))F{A = 0) , 
\\H{si) - H{s) - (si - s)H{s)\\ < C\si - sf , 

see (jl.4j) . Assumption Al and Uq = F{A = 0)IIq, we are left with showing that 

f dsi\si-s\'^-\\F{A = 0)us{si-s)I±{s,a)\\i<Ce . (4.9) 



We insert 1 = F{A < a + a) + F{A > a + a) to the left of the commutator in (jl.lH) . By 
(|3.3l 13. 8|) the trace norm appearing under the integral is bounded as 

||F(A = 0)u,(si-s)/±(s,a)||i 

< C\\F{A = 0)us{si - s)xiH{s))F{A < a + a)|| + Cn{1 + a)"^ , 

where we also used ()3.6p . By ()3.2j) with b = the latter norm is estimated as 

\\F{A = 0)us{si - s)x{H{s))F{A < a + a)\\ 

< C\\F{A = 0)F{A <a + a + ee-\si - s))\\ + Cn{1 + - s])"^ . (4.10) 
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Picking a = \9e ^{s — Si) > 0, the first term vanishes for s — Si > 2a6 ^e. Now Eq. (j4.9|) 



2 

holds because the I.h.s. is estimated by a constant times 



_-2 



poo 

' / da a\F{a < 2ae-^e) + (1 + V)"^) < C{ea^ + e) , (4.11) 
Jo 

which proves the lemma. □ 
For later purposes we retain the following estimate from the above proof: 

\\FiA = 0)e-'^(^)*/±(s, a)|| < < 2ae-') + C^(l + . (4.12) 

In particular, it proves the first remark after Prop. 14.11 

The last term in Eq. ()4.8|1 is also 0e, as shown by the next estimate. 
Lemma 4.3. We have 

\\X,{s.,s)I±{s,a)X,{s.,sy\\i < Ce^ . (4.13) 
Proof. As in the proof of the previous lemma, the norm in ()4.13p is bounded by 



s 



dsids2\si - s\\s2 - s\\\F{A = 0)us{si - s)I±{s,a)usis2 - syF{A = 0)||i . (4.14) 

Using again (j3.3l 13. 8p . this last norm by itself is bounded by 

(\\F{A = 0)us{si - s)x{H{s))F{A < a + a)|| + Cn{1 + a)""^) • (si ^2) 

where (si — > S2) is shorthand for the previous expression with si replaced by S2- We pick 
a = \0£~^ min(s — Si,s — S2) > 0. Then the previous expression is estimated as 

(^F((Ti < 2ae-^e) + (1 + ^Oe-^ min(ai, (72))"^) ■ (si S2) , 

where (Tj = s — Sj. Using 

f°° 1 

/ da2 ^2(1 + -ee-^ min((Ti, ^2))"^ < C{al + e^) , 
Jo ^ 

/ daido2 (Ji02{l + 7.(^^~^ min(cri, a2)Y^ < Ce^ , 
Jo ^ 

the claimed bound is established for ()4.14|) . □ 

We now turn to the second and third term in Eq. ()4.8|) . which will account for the 
current in the adiabatic limit. 

Lemma 4.4. 

lim II (f/,(s_, si) - Usis. - si))usisi - s)piH{s))xiH{s)){±A - a + = , (4.15) 

uniformly also in si G [s-, s]. 
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Proof. We first consider the case where p is smooth. Since T = x(if(s))(±A — a + i) ^ G 
Ji by (TOH]) and s-lim„^oo ^(^ > «) = we have hm^^oo > a)T\\i = 0. By 

approximation we may thus prove 

hm \mis., si) - Usi-s^ - s,))usis, - s)piH{s))xiH{s))F{A < a)|| = , 

where x ^ C*^(0, oo) with xx = X- To estimate this operator norm we use ()4.6|1 together 
with 

\\H{s2)-H{s)\\<C\s2-s\ (4.16) 

to obtain the bound 

e-' ^ ' ds2\s2 - s\\\F{A = 0)m,(s2 - s)p{H{s))x{H{s))F{A < a)\\ < Ce{l + , 

where we used, see the argument in connection with ()4.10j) . that the norm under the 
integral is bounded by F(s — S2 < e9~^a) + CN{i + £~^\s2 — s\)~^. The general case, where 
p is of bounded variation, is also dealt with by approximating p{H) = s— lim„^oo Pn(-f^) 
with pn smooth. In fact we can pick n so that — p{H))T\\i is arbitrarily small 

(uniformly in s). □ 

Proof of Prop. I4!T1 As mentioned in the introduction, [/^(s, s_)p(iJ_)t4(s_, s) is inde- 
pendent of S- < 0. We may thus evaluate the r.h.s. of ()4.8|1 in the limit s_ —>■ oo. By 
estimate and its adjoint, this amounts to replacing p{HJ) by p{H{s)), i.e., 

e-' tT{p{H.)Ue{s., s)Ii{s, a)Us{s, s.)) = e'' tT{p{H{s))I±{s, a)) + 
Lim tj:{Xe{s^, s)*Us{s^ — s)p{H{s))I±{s, a)) + 

s-^— oo 

Lim tr(u,(s - s_)Xs(s_, s)/±(s, a)p(/7(s))) + 0e , (4.17) 

where we also used [us{s — S-),p{H{s))] = and ()4.13p . The first term one the r.h.s. is 
just the equilibrium value of the current subtracted on the l.h.s. of Eq. ()4.2|) . As we are 
not going to show that the limits in ()4.17p exist, they are just to be understood as sets of 
limit points. We next claim that 

e-^ Lim tT{{X,{s-,syUs{s- - s) -en^l\s))p{H{s))I±{s, a)) ^ , {e ^ 0) . (4.18) 

By taking the complex conjugate this implies 

e-^ Urn tr{I±{s,a)p{H{s)){us{s - s^)Xeis^,s) + en^l\s))) ^ , {e 0) , 

s_— >oo 

where we used fi^''(s)* = — i7^''(s). Used together in ()4.17|) they prove ()4.2|) . 
We next note that, by a change of variables, 

efi^Vs) = —ie~^ / dsi(si — s)us(s — si)H(s)us(si — s) . 
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As noted before, the integral is convergent in trace class norm upon multiplication from 
either side with I±{s,a). For the purpose of proving ()4.18p we may thus assume it to 
have s_ as its lower limit of integration. Then the trace there, including the factor in 
front, equals 



— le ^ J dsi{si — s)x 

tr{us{s - Si)H{s) (f4(si, s^)us{s_ - s) - Us{si - s))p{H{s))I±{s, a)) . 

We use 

UeiSi,S^)Us{s^ - s) - Us{Si - s) = [/^(Si, S_) - Si) - [/^(s. , Si)) - s) 

and turn Us{s — si)H{s) around the trace, so that the previous expression estimated as 

e^^J dsi\si - s\\\{Ue{s^, si) - Us{s^ - si))us{si - s)p{H{s))x{H{s)){±A-a + iy^\\iX 

\\{±A - a + if[H,{s), f{±A - a)]x{H{s))us{s - s^)F{A = 0)|| . 
The first factor tends to zero uniformly as e ^ by (j4.15p . so we are left to show 



J dsi\s, - s\\\F{A = 0)us{s, - s)x{H{s))[Hk{s), f{±A - a)]{±A -a + i)^\\<C . 

(4.19) 

We insert 1 = F{A <a + a) + F{A > a + a) to the left of the commutator. By (jSSl 
and the remark following them we obtain the bound 

£-2^ dsi\si - s\{\\F{A = 0)us{si - s)x{H{s))F{A < a + a)\\ + Cn{1 + a)-"") , 

where we take a = ^9e~^{s — Si) > as in (j4.1Up . The resulting bound differs from (j4.11|) 
by having a instead of cr^. Hence the bound (j4.19p . □ 

5 The scattering limit 

In the previous section we saw that the currents can be computed from frozen data in the 
adiabatic limit. These data were not directly related to the frozen scattering data and 
involved objects like H{s) and fi^''(s). In this section we show that in the limit that the 
ammeter is far, a —>■ oo, then all we need to know is the frozen, scattering operator S and 
the initial state p. 

First we show that in the large a limit the incoming, resp. outgoing currents have no 
scattering in the past, resp. the future. As usual, all statements are uniform in s G /. 



Lemma 5.1. We have 



lim \\n^^\s)I^{s,a)\\i = , 

a^oo 

hm i|(fi«(s)-S«(s))/+(s,a)||i = 0. 
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We recall that S^^^ was defined in (j4.5p . These facts will yield a first expression for 
the scattering limit of the current, Eq. ()4.4j) . 

Lemma 5.2. We have 

hm tr([fi«(s), p(/7(s))](/+(s, a) + /„(s, a))) = tr ([^^(s), p(iJ(s))]/+(s, a)) , (5.1) 

a— >oo 

where a on the r.h.s. is arbitrary. 

We then express the latter result in terms of "frozen" scattering data, such as the 
scattering operator S{s), defined in Notice that it acts on the asymptotic Hilbert 

space L^(R+,C") of the channels, rather than to H. Further distinguished operators 
acting there are the Neumann Hamiltonian Hq, introduced in the introduction, and the 
generator of dilations, 

I d d 
2i dx dx 

which may be regarded as a model for ()1.6|) . The trace ()5.H) can then finally be computed 
exactly. 

Proposition 5.3. Suppose, in addition to the hypotheses of Thm. \l.ll that p is smooth. 
Then, 

tr{[S^'\s),piH{s))]l4s,a)) = -itr(5(s)5(s)V(/7o)n,i[ifo, /(^ - a)]) (5.2) 



(5.3) 



where Aqj = AoUj, and the scattering operator S{s) is defined in Eq. / Ii..9)) . 

For p which is not smooth an approximation argument will complete the task. 

Proof of Lemma f5.1[ Both claims of this proposition are the consequence of the stronger 
bound: 

||n£Hs)/±(s,a)||i<C^a-(^-2), 

for large enough A^. Let us proof this bound for a case of the proof for follows 
the same lines. It is clear from Eq. ()4.3j) that we may establish that bound for 

/ dt\t\ ■ ||F(A = 0)e-'^(^)*J_(s,a)||i . 

J — oo 

The norm under this integral also appeared under the integral (j4.9p . We estimate 
it similarly by means of (j3.3l 13. 8|) . except that we insert F{A < —a + a) (and the 
complementary projection) instead of F{A < a + a). By ()3.2j) we get 

||F(A = 0)e"'^(")*/_(s,a)||i < C\\F{A = 0)F{A < -a + a - b + et)\\ 

+C^(|6| + |t|)-^ + Cjva-^. 

Choosing b = —a/2, a = (a/2 — 6t)/2 > 0, we see that the first term vanishes and we 
obtain the desired bound since 

fO 



[ dt\t\{^ + e\t\)-'' <CNa- 

J -oo 2 



(7V-2) _ 

□ 
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An ingredient to the proof of Prop. 15.31 and, to a minor extent, of Lemma [5.21 are the 
relations Eqs. (ITTKl ^J^. Formally, Eq. (ITTHll follows for p(A) = e^^*, (t G R), by a 
change of variables in ()1.14p and hence for general functions p. For our purposes we shall 
need a somewhat more precise statement, given by the first part of the following lemma 
(cf. the first remark after Prop. I^TTjl . 

Lemma 5.4. (a) Eq. M.l(^) and, for p smooth, Eq. / li. j,5|) hold true if multiplied from 
either side by a). 

(h) The wave operators Vl±{s',s) = ?,—\im.T^±oo^^^'^^''''^^~^^'^^'''^Pa.c{H{s)) for the pair 
{H{s'),H{s)) satisfy the intertwining property 

f{H{s'))n^{s', s) = n^{s', s)f{H{s)) (5.4) 

for any (Borel) function f . 

(c) The scattering operators for the pairs {H{s), H{si)), {i = 1,2) are related as 

S{s, S2) = ^+{s2, si)S{s, si)Q^{s2, si)* . (5.5) 

Proof. Parts (b, c) are standard |221 ISj- The wave operators exist and are complete 
under our assumptions. The chain rule for wave operators, 

Q±{S, Si) = Q±{S, S2)^±{S2, Si) , 

and S{s,Si) = fl+{s, Si)*fl-{s, Si) imply S{s,si) = Q-^{s2,si)*S{s,S2)fi-{s2,si). From 
this ()5.5|) follows by the completeness of scattering, Q±{s2, Si)*fl±{s2, Si) = Pac(-f^("S2))- 

We begin the proof of part (a) by claiming that 

\m4s\s)-l)x{H{s))\\<C\s'-s\, (5.6) 

where x e C^(0,oo). Indeed, let x e C^(0,oo) with XX = X- Then, by (Q, - 
l)x{Hism < \\xiHis'))in+ - l)xiHi.sm + \\xiHis'))xiHi.s)) - xiHis))]]. The second 
term is bounded by \\x{H{s')) — x{H{s))\\, and fits the bound ()5.(')|1 by ()A.5j) . For the first 
term we use the fundamental theorem of calculus: 

x{His'm^is',s)-l)x{H{s)) = s-]im 1 C dtx{H{s'))e^''^^'^'Ve-^''^^^'x{H{s)) , (5.7) 

where V = H{s') - H{s). We write V = {AyViAy and use KU^ for V, as well the 
if(s)-smoothness of {A)~'^x{H{s)) for r > 1/2 (and similarly for s'), see item (ii) in Sect.|21 
An application of the Cauchy-Schwartz inequality on matrix elements of ()5.7p yields ()5.6p . 

We can now prove the statement about Eq. (jl.lfij) . whose r.h.s. we denote by T. Then 
TJ+(s, a) is a convergent integral by ()4.12|) . Moreover, Te~^^^'^'^^I^{s, a) = e~^^^'^^^TI^(s, a), 
whence TxiH{s))I+{s, a) = xiH{s))TI+{s, a) for x e 0^(0, 00). Since xiH{s))I+{s, a) = 
/_i_(s, a) for XX = X and x as in ()1.11|) . we may thus prove 

/oo 
c?te^^(^)*if(s)e-^W*J+(s,a) . (5.8) 
-00 
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We write S - 1 = - n+), so that 

/oo 
rfte'^(^')Ve-'^W*J+(s,a) 
-oo 

/oo 
dte'''^'^'xiH{s)){n\ - l)l^e-'^(^)*/+(s,a) 
-oo 

/oo 
rfte'^(^)Ve-^^(^)*/+(s,a) . 
-oo 

The first integral is estimated as C(a)|s' — due to Eqs. ()5.6l I4.1bl 14. 12^ . In the second 
one, the contribution coming from the remainder in = H{s){s' — s) + 0{{s' — s)^) is 
estimated the same way. This imphes ()5.8|1 . 

As for Eq. ()1.15p in the case p{X) = e^'^*, the change of variables mentioned before 
is legitimate, as the integrals are convergent once multiplied with J+(s,a) due to ()4.12p . 
The result extends to p{H{s)) with p G C^(M) or p' G C^(M), which amounts to the 
same since H{s) is bounded below. □ 

Proof of Lemma f5.2[ By Lemma I^TI the l.h.s. of ()5.1|) equals lima^oo ^^{[S^^\s), p{H{s))] 
/+(s, a)), provided this limit exists. It does, since the expression is independent of a. For 
p smooth this follows from Lemma 15.41 Indeed, the r.h.s. of ()1.15j) commutes with 
H{s), so that the independence is seen as in ()2.H|1 . For general p the conclusion follows 
by approximation by a sequence {p„} of smooth functions, such that s— limp„(if(s)) = 
p{H{s)). The traces converge by ()3.19p . □ 

Proof of Prop. 15. HI For the sake of precision we recall that the scattering operator S{s) = 
depends on wave operators defined in a two-Hilbert space setting: n±(s) = 
s-limT^±ooe'^(')^Je-'^o'^, where J : L'^{R+,C'^) ^ H is the embedding given by (IT^ . 
We may also replace J by a smooth function J on R, which equals near a; = 0, and 1 
near x = oo. This is without effect on r2±(s), since s—limT^±oo{J — J)e~^^'^^ = 0. In 
particular, we may pick J so that J = 1 on suppv, see In the present context (|5.5p 

reads S{s',s) = Q+{s)S{s')Q^{s)* and implies 

ds,s{s',s)s'=s = ^-,{s)s{s)n4sy , 

where the differentiability is granted e.g. after multiplication with /+(s,a). 

We now prove Eq. ()5.2|) with Ilji[Ho, /{Aq — a)] replaced by Ilji[HQ, f{A — a)] = 
i[Ho, f{Aj — a)]. Its l.h.s. equals, by Lemma EiH 

tT{n4s)S{s)n4srp'{H{s))x{H{s))l4a)x{H{s))) 

= tT{S{s)S{srp'iHo)x{Ho)lMxiHo)) 

+ tT{S{s)n4srp'{H{s))x{H{s))[l4a),n+{s)]x{Ho)) . 

Here /+(a) is defined in Eq. ()1.7|) . We turned around the trace, which thereby 

moved from Ti to L^(M_|_,C"). The first term corresponds to the claim, since xp'x = p'- 
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It is again independent of a. The second term will thus vanish as soon as it does in limit 
a oo. The commutator there is [/+(a), = [/+(a.), — J] by the above choice 

J. We are so left to show that 

hm \\x{H{s)){n+is) - j)xmiM\u = , 

a— >oo 

and also with fi+(s) replaced by fi+(s)* and x{Ho), x{H{s)) interchanged. Using the 
integral representation for — J the estimate reduces to 

POO 

/ dt \\xiHis)[Ho, J]x(/^o)e-^^«*/+(a)||i 
Jo 

POD 

<C dt \\FiA = 0)x(i/o)e-'^»*/+(a)||i ^ , (a ^ oo) , 

^0 

(and with Hq, H{s) interchanged) which holds true by the estimates in the proof of 
Lemma 15.11 

We next replace A by Aq on the r.h.s. of ()5.2|1 . Since both traces are independent of a 
we may, in each of them, replace / by a smeared switch function /', such that / G H'^{d) 
with d> 2 (see the Appendix for notation). Since {f{Aj)-f{Aoj))p'{Ho) e Ji by (IA.12jl . 
the difference is seen to vanish by cyclicity. 

Finally, to prove ()5.3j) . we may keep / replaced by / as before. It follows from 
Eqs. (|Xl0l 1X21) that 

tT{S{s)Sisyp'{Ho)n,[Ho, f{Ao - a)]) 

= tT{S{s)S{srp'{Ho)U^{f{Ao - a - 2i) - /(Ao - a))Ho) 

1 roo 1 rp / fl ^ \ 

To compute the last integral, note that f^^dX (/(A — c) — /(A)) = — c for any c G R 
which extends by analyticity to | Imc| < d. Hence the result. □ 

Proof of Theorem ll.il We have shown that 

tT{[S^'\s),p{H{s))]l4s,a)) = -^l°^ME){^S*)^^ (5.9) 

holds true if p is smooth. If p is of bounded variation it can be approximated by a 
sequence p„ G C^(]R) with Pn(A) p(A) pointwise and uniformly bounded total variation 
sup„V"(p„) < oo. Then s— lim„ p„(if(s)) = p{H{s)) and dpn{X) — > dp{X) in the sense of 
weak* convergence. Hence ()5.9|) is inherited by the limit. The proof is completed by 
Eqs. dOI) and □ 
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A Appendix. Trace class properties of g{H)f{A) 



In this section we will discuss different properties of the operator product g{H{s))f{A). 
We first prove that operators of the form g{Ho)f{Ao) are Hilbert-Schmidt under suitable 

conditions on the functions g and /. For operator of this type one can essentially write 
down the integral kernel, and use it in order to compute the Hilbert-Schmidt norm. 
Heuristically the trace of an operator is the integral over phase space of the symbol 
(divided by 2%). Since 

under the map {x,p) ^ {E — p'^,a — px), we introduce the norms 

/oo poo r r-i 

Ja\fia)\\ \\\g\r^ ^1^(^)1'' 

on functions / on the real line, resp. g on the half-line. The Fourier transform is f{t) — 
(27r)-Ve-'^7(A)dA. 

Proposition A.l. We have 

UH,)f{A,)\\, = {2n)-"^\g\\\-\\f\\ . (A.l) 
If g and f are continuous and the operator gf(i7o)/(^o) ^■s trace class, then 

tr {g{Ho)f{Ao)) = -J^ ^9iE)-J ^ da f{a) . (A.2) 

Proof. The space L^(]R+) can be identified with the even functions in L^(R) by means 
of the map {Jip){x) — (p{\x\), with || Jv^|P = 2||(/9|p. Consider the operator g{p'^)f{Ao) on 
L^(R). Since it preserves parity and reduces to Hq on the (so identified) subspace of 
even functions (with projection P+) we have 

\\g{Ho)fiAo)\\l = tr^.(K^)(^(//o)|/r(^)^(/^o)) 

= ^tr^.(«)(P+^(p2)|/|^(Ao)<7(/)P+) 
= ltrr^^u)iU*P+g{p')mAo)gip')P+U) 

= ltrr..^^){P+g{x')\m-Ao)g{x')P+) 

= trr^.^^^){-g{x')m-Ao)g{x')) , (A.3) 

where {Ui)){k) = (27r)i/2^(A;) is the Fourier transform: [U,P+] = 0, P+ = (1 + U^)/2, 
pU = -Ux, and AqU = -AqU. Using the MeUin transform M : L'^{R+) ^^(R), 

(M<^)(a) = (27r)-V2 / -—x'-<f{x) , 
Jo ^ 
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which diagonahzes Aq = M*aM, one obtains that h{a) has integral kernel 

hiA)ix,y) = (27r)-^(xy)-U(log^) (A.4) 



X 



for h e L^(R). Since the kernel is continuous in x,y > 0, we can write for ()A.3jl (cf. 
Thm. 3.9. [211) 

i2n)-' r-\9ix')\' r dam-a). 

Jo ^ J~oo 

By means of the change of variable = E, this is the square of the r.h.s. of (|A.1|) . 

Using the Fourier transform as in ()A.3|1 we find tr((7(ifo)/(^o)) = tr((yf(a;^)/(— Aq)) . 
The integral kernel of the latter operator is obtained from ()A.4j) and the trace is the 
integral over its diagonal, due to Thm. 3.9 ji24j. □ 



In the rest of this appendix we shall give an example of trace class operator involving 
Aq, Hq and use it to estabhsh g{H{s))f{A) G jTi for a large enough class of operators. 

Lemma A. 2. Let g e C^O^), then 

\\g{H{s))-g{H{s'))h < C\s - s'\ , (A.5) 
g{H{s))-JgiHo)r G Ji , (A.6) 

where J : L^(]R+,C") Ti is the embedding given by il.^) . 

Proof. We first prove ()A.5jl for g{X) = (A + i)"^*", where m is as in Assumption A2. 
Setting R{s) = {H{s) + we have by ()1.4j) and the resolvent identity 

2m 



fc=l 

2m 



k=l 



The desired bound now follows from Al by using A2 either for s or s'. For general 
g G C^(M), as well as for the rest of this proof, we use the Helffer-Sjostrand representation, 

g{H) = - [ {H - z)-'d,~g{z)dxdy , (A.7) 

where = [d^ + idy)/2, and g is an almost analytic extension of g with g{z) vanishing to 
a high power near the real axis. Before doing that, we set G(A) = 5'(A)(A + i)^™", so that 
the first term in 

giHis)) - g{H{s')) = - i?(.')'")G'(i/(.0) + - 

is taken care of. For the second, we apply ()A.7|) to G and are led to estimate 

\\R{sr"^[{H{s)-z)-'-{H{s')-z)-']\\,< 

\\{His) - z)-'\\\\Risr"^UMHis') - His)\\UH{s') - z)-'\\ < C\ Imz|-2|s - .'| , 
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which completes the proof of ()A.5|) . 



Before proving ()A.6|) we claim that 

9{H)h G Ji , (A.8) 

where h acts as multiplication by /i G C^[0, oo) on L^(M_|_,C"') and by /i(0) on TYq- 
Since g{H)IlQ G Ji it will be enough to show (1 — YlQ)g{H){\ — Yiojh G Ji or actually just 
\\J{H — z)~^ J*h{x)\\i < C|Im2;|~^. Since the kernel of {H — z)~^ is decaying exponentially 
0, matters are further reduced to 

\\Jl{H - z)-^Jl\\i < C\lmz\-^ , 

as an operator on L^([0, L], C"), where Jl is the corresponding embedding operator into 
Ti. The initial piece [0, L] has to be taken large enough, so that supp h C [0, L). Let now 
B be the quadratic form of the Bilaplacian, 

B{ip,iP)= [ ^"{x)i)"{x)dx , (A.9) 

with domain given by the Sobolev space W^^(0, L). We maintain that 

11(1 + Bf/^JL{H - z)-^Jl\\ < C\ lmz\-^ , 
(1 + B)-^l^ G Ji , 

which imply the claim. The first statement follows by A3 and ()A.9|) from 

Jl{H - z)-^JIBJl{H - z)-^Jl = T*T , 

with T = 1 + zJl{H — z)'^ JI- The second statement can be seen by computing the eigen- 
values /c^ of the operator B associated to (jA.9jl . The latter is given as -B = d^/dx^ with 
boundary conditions if" = (p'" = at x = 0, L. From this one computes the eigenvalues 
as the zeros of 1 — cosh kL cos kL. 

To prove ()A.6|) we make use of a smooth embedding J as in the proof of Prop. 15.31 
Since (J — J)g{HQ) G J^i (see j2I], Thm. XI.20) we will prove the trace class property for 

g{H) - Jg{Ho)r = g{H){l - J J*) + {g{H)J - Jg{H,))r . 

The first term is trace class by ()A.8|1 . For the second term, Eq. (1A.7j) and the resolvent 
identity reduce matters to 

\\{H - z)-'[J,H]{Ho - z)-'J*\\, <\\{H- z)-'[J,H]\\\\hix)iHo - z)-'\\, < C\ lmz\-' , 

where h G C^{R) with hJ' = J'. □ 

The Hardy class H'^{d) consists of all functions /, analytic in the strip {z \ \ lmz\ < d}, 
which satisfy 

/oo 
\f{x + iy)\'^dx < oo . 
'OO 

We recall that / G H'^{d) is also characterized by supo<y<(i ||e^''^'/(^)||2 < oo. 
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Proposition A. 3 (The pull-through formula). For any f E H'^{d) with d > 2 

HofiAo) = f{Ao - 2i)Ho . (A.IO) 

Proof. We first prove the statement for f{x) = e~'*^'. Under conjugation with the dilation 
operator, (e~^^'^°tjj){x) = e~^^'^xjj{e~^x), the Neumann Laplacian becomes 

e'*^»i/oe"'*^° = e-^'Ho , (A. 11) 

hence i7oe"'*^° = e-'^^^^'^'^Ho. Now, for / G H'^id), 

/oo poo 
H,e^'^^f{t)dt = / e'*(^°-2')i7o/(t)c^t = /(Ao - 2i)Ho , 
-oo J — oo 

where the last step is justified by the above mentioned property of the Hardy class func- 
tions. □ 

Lemma A.4. Let f G H^{d) with d > 2 and g e C^(M). Then 

{f{A)-f{Ao))9{Ho)eJ,. (A12) 
Proof. As mentioned, the Fourier transform of /'(A) = Jf^ dte^^^f'{t) is bounded as 

\f'{t)\<Ce-y\'\ , (A13) 

for any y < d. We represent / as 

/•oo .QiXt_i.^ 

/(A) = /(o) + j jt[^-^)m, 

and write 
1 



Jo 

which follows from ()A.11|) . Applying Thm. X1.20 21J to 

{A - Ao)g{Hoe-'') = {w{x)p - '-w'{x))g{H,e-^^) , 

where w{x) = v{x) — x, one finds 

\\{A-Ao)g{Hoe-'^)h < Cs{e^il+^^ + e^^'^+'^) , (A.14) 

for any 5 > i. We obtain \\t-\&^^ - e'^«*)^(i/o)||i < C{1 + e^^) with d= l + S. Picking 
6 so that d < d, we obtain the claim from Eq. ()A.13|) . □ 

The following example will be useful. 

Lemma A. 5. 

\\Ho{Ho + l)-^{Ao - z)-^\\, < C\ lmz\-^ . (A.15) 
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Proof. It suffices to prove the claim for | lmz\ large enough since \\{Aq — z + iy){AQ — 
^ C\ Imz|~^ for small | Imz|. We shall do that for | Im2;| > 4. By the pull-through 
formula, Eq. ()A.10|) . we compute 

Ho{Ao-z)-' = {Ao-z -2i)-' Ho 

= {Ao-z- 2i)-\Hl + 2Ho + l)Ho{Ho + 1)"' 

= [Hl{Ao -z + 2i)-i + 2Ho{Ao - z)-^ + {A^ - z - 2i)-']Ho{Ho + l)"' . 

We multiply this expression from the left by {Hq + 1)"^ and from the right by {Aq — z)^'^, 
so as to obtain the bound 

1 

\\Ho{Ho + l)-\Ao - z)-% <J2\\{Ao-z + 2ki)-'Ho{Ho + l)-^^ - z^rX . 

k=-l 

Eq. (jSm now yields \\hI^^{Ho + iy\Ao - z)-^\\2 = (2| Imz|)-i/2_ YLence f)Xl5|l follows 
from the Holder inequality. □ 

Lemma A. 6. For g G C^(0, oo) we have 

snp\\{A-a±i)-^g{H)\\i<oo . (A.16) 

In particular, for any f G C°°(R) with f{x){x + i)^ bounded, 

f{A)g{H) G Ji . (A.17) 

Proof. Since g{HQ)HQ^{HQ + 1)^ is bounded, by ()A.15|) the bound holds for Aq, Hq in 
place of A, H. We ffist undo the replacement in A and write 

(A - a ± ir'g{Ho) = {A-a± ir^<P{x)g{Ho) + {A - a ± - m)9iHo) , 

where G C^[0, oo) has v{x) = x on supp(l — 0). Then (j){x)g{HQ) G Ji (see [21], Thm. 
XI.20). For the second term on the r.h.s. we write, dropping a, [A + — (j){x)){AQ + 

i)^ ■ {Aq + i)~^5f(iJo)- Since the last factor is trace class, we are left with showing that the 
ffist factor is bounded. This follows from 

(1 - 0)(Ao + = iA + 1)2(1 - 0) - (A + i)[0, Ao] - [0, Ao]iA + i) , 

where [0, A] = ix(f)'{x) bounded. 

The bound (IXl6ll clearly also holds for {A + i)-^Jg{Ho)J* = J{A + i)-^g{Ho)J* and 
hence, by ()A.6|) . as claimed. □ 
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